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Abstract 

We develop a method for constructing asymptotic solutions of fi- 
nite - difference equations and implement it to a relativistic Schroedinger 
equation which describes motion of a selfgravitating spherically sym- 
metric dust shell. Exact mass spectrum of black hole formed due to 
the collapse of the shell is determined from the analysis of asymptotic 
solutions of the equation. 



1 Relativistic Schroedinger equation in finite 
differences 

In this note we develop a method for determining the energy spectrum of 
relativistic Schroedinger equation in finite differences 

*(m, fi,S + i() + *(m, /i, 5* - z£) = 

(F m F out y 1/2 (F m + F ou4 - M 2 /4m 2 5) tf (m, jx, 5) (1) 

which appear to describe a quantum black hole model with selfgravitating 
dust shell introduced in [|1[] . Here parameter M is the bare mass of the shell, 
Tn = m out is the Schwarzchild mass of the space-time outside the collapsing 
shell, /j, = m in /m out is the ratio of Schwarzchild masses inside and outside 
the shell, S = R 2 /4G 2 m 2 (G is gravitational constant) is the dimensionless 
area of the shell surface and F inout = 1 — 2Gm in ^ out / R are functions entering 
the Schwarzchild line element. 

Parameter ( = m p (/ (2m 2 ) where m p i = J chj G is Planck mass (h is 



Planck constant and c is the speed of light). It becomes small in semiclassical 
(H — > 0) limit. In this limit the finite shift in the argument of the wave 
function in ([]]) become small and one could us an approximate expression 

^(S ± <) « *(£) ± i(*'(S) - ^*"(S) (2) 

Then the equation [I] becomes just a usual Schroedinger of nonrelativistic 
quantum mechanics: 

- C V + {2 - (F m F out y 1/2 (F m + F out - M 2 /Am 2 s) } * = (3) 

This equation was studied in |l|, in details. There was an expression for 
the energy spectrum of bound states found and quasiclassical behaviour of 
the wave function studied. 

In this paper we will analyse the solutions of the exact equation (|]). 

The shift of the argument in the wave function is along the imaginary 
line, so the equation is naturally defined over some complex one dimensional 
manifold. The function 



(F- F .W 2 

V- 1 in ± O' 



out) 



,,,--, igW'"- 1 "'-") (4) 
\ V pJ V pJ p 




(p = yS) which enter the equation ([!]) is a branching function of its argu- 
ment. So ([J) is defined over the Riemannian surface Sf for the function (^) 
which is a two dimensional sphere glued from the two Riemannian spheres 
S + and S- along the sides of the cut made along the interval (//, 1) of the 
real line. 

There are two points p = oo on the Riemannian surface Sf in S+ and S- 
components, which are singular points of equation ([!]). They differ by the 
different chooses of the sign of the function {F in F out ) 1 / 2 . 

The natural requirements on the wave function \l/ of a bound state is that 
it must decrease with p — > oo But in our case we have two infinities in S + 
and S- components of Sf- As it was argued in [p], |, §] the two possibilities 
for p — > oo have the following physical meaning. If one consider maximally 
extended spherically symmetric solution of Einstein equations in vacuum - 
Kruskal manifold M, one finds that the radial coordinate p take the same 
values on different sides of Einstein-Rosen bridge (in R+ and i?_ regions 
of Kruskal manifold [0). So p could tend to infinity twofold: either in R + 
region or in i?_ region. This complicated structure of Kruskal manifold is 
mirrored in complicated structure of configuration space of quantum system 
(the real section of Sf) which describe a self gravitating dust shell. So the 
two possibilities for p — > oo in S+ and SL components of Sf corresponds to 
the possibility for p — > oo in R + and R- regions of Kruskal space-time. 

But the requirements for the wave function to vanish at infinities are not 
all the requirements which a reasonable wave function of a physical state 
must satisfy. Besides if must be one- valued function. 

This requirement is satisfied automatically if we consider Schroedinger 
equation over a complex plane (or over the real line). But the equation ([I]) is 
defined over (the real section of) the Riemannian surface Sf- The two infini- 
ties in S± components are both singular points of equation. So the equation 
([!]) is defined in fact over complex manifold which is two-dimensional sphere 
minus two points - a cylinder. But a cylinder possess a nontrivial cycle. The 
monodromy of solution along this nontrivial cycle is not identity in general. 
Therefore, the requirement for the wave function to be one- valued function on 
the Riemannian surface Sf (and on its real section) become a nontrivial one. 
In what follows we will find the asymptotics of exact solutions of equation (|l|) 
at both infinities. The analysis of analytical properties of asymptotics and 
the requirement for one-valued solutions to decrease at infinities will lead us 
to defining of the mass spectrum for the equation ([!]) 



2 Asymptotic solutions 

Let us find the asymptotic solutions of this equation at infinities in S± 

components of Sp- 
it is convenient to rewrite ([3]) in matrix form. 
If we truncate the Tailor expansion 

27V (_Af\k 

*(s - K) = * + £ Hr"* w ( 5 ) 

fe=i ft! 

on 2iV-th item and expand asymptotically 

(-*r(-Ar- 

1 + /i 1 3/x 2 + 2/i + 3 1 5/i 3 + 3/i 2 + 3/i + 5 1 
+ ^~7S + 8 S + 16 ^ () 

in powers of S, we could rewrite ([[J) as 

«) 2 ^ 2- (2k)\ 

(gjO! f_i + 1 + (W) 2 -^ 2 /^ 2 + 
«) 27V V 85 

(l + / i)(2(l-/i) 2 -M 2 /m 2 ) \ 

16s 3/2 J* IM 

where the upper signs correspond to 5 + and the lower ones to £_ components. 
It is convenient to introduce a vector function 

y 1 = tf; Y 2 = V'; ...; Y 2N = V™' 1 (8) 

Then (0) takes the form 

Y' = AY (9) 

where 

A(S) =A + -LA! + ±A 2 + ^A 3 (10) 



where Aq has the following form: 



A 



( ° x 






1 


. 
. 






0\ 





^ oi 



a 3 


. 




■ 02JV-1 


1 



with coefficients 



(2iV)! , 



02fc+l 



(2JV)! 



Matrices ^2,3 are 
and 



(2fc!)«) 2JV - 2fc ' 



a 2 fc = 
£;>0 



A 2 = a 



^ ... ^ 
... 

V i o ... o j 



a 



{2N)\ {{I - pf - M 2 / m 2 ) 
«) 27V 8 



/ ... \ 



(11) 



(12) 
(13) 

(14) 



A, = (3 







{2N)\ (1 + n) (2(1 - /i) 2 - M 2 / m 2 ) 



«) 



2 A' 
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... 

Vi o ... oy 

After the change of variable S = p 2 equation @ takes the form 

-Y; = AY; A(p) =A + - 2 A 2 + -~A Z 
p r p z /r 



(15) 



with matrices 



Ai — 2iAi 



(16) 



(17) 



If we let N — + oo then matrix equation (|9|) is equivalent to finite - differ- 
ence equation ([[]) in asymptotic regions. There exist well-developed method 
for finding asymptotic solutions (p — > oo) of matrix equations of (0) type [H]. 



First we need to make a transformation 

Y = TY 
A = TAT' 1 

Y p = ~AY (18) 

such that matrix A becomes diagonal: 

A = diag(\ 1 ,...,\ 2N ) (19) 

The eigen values Ai,...,A 2 at could be found by equating the characteristic 
polynomial for matrix Aq to zero: 



P2N 
As a result we get an equation 



A n -XI 



(20) 



^-^(i^(?r ± fo (21) 

In the limit iV — > 00 this equation is equivalent to the following one 

cos ( ^ = ±1 (22) 

So we see that all the eigen values become double degenerate in the limit 
iV —>■ 00 as it is seen from Fig. |l|. This means that the matrix An could not 
be transformed to diagonal form in this limit, instead it could be transformed 
to Jordanian form: 

An = TAnT- 1 = ,7^.12® ... © J N ; (23) 

Ji = AjJ + H\ 

/= (o i) ; ^ = (o 0) 

It is convenient to rewrite all the matrices in the form of blocks of dimension 
2x2: 

r =C -~ T J\ r -=(::;:;:) <*> 

J-Nl ••• J-NN J 



COS 



(?) 



Figure 1: Double degeneracy of eigen values mJV^co limit 



and in analogous way 



A 



u 



2 




A 



NK 




2a 2 K-i 



X 




(25) 



Then the matrix T could be found from the equation ( |23| ) which lead to the 
following system of recurrent equations: 



Solving this system of equations gives T in the following form 



(26) 



T, 



i.i 



(¥) -%w (v)"-*4 J, + </-i)(¥)' M 4M 



/ f^ 



U ; 



ii 
l 



n 



27- 1 



(j) i 2/-i r^ 2(J ~ i:) ^ j ) 



t 2 + 



(27) 



5ii' / 



the coefficients 4i an d A2 are arbitrary and the form of asymptotic solution 
of equation (^) do not depend on them. So we will take % = t^ 2 = 1 f° r 
simplicity. 
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Matrices A 2 a transform under CIT8D as follows: 



A, 




(28) 



where t\- are elements of the block (T r ) IN of inverse matrix of transforma- 
tion T- 1 . 

The transformed matrix equation (p~8| ) could be reduced to a set of 2 x 2 
matrix equations with the help of one more transformation 



where 



Y = P{p)Z 



( P 12 

P21 



P2N 



(29) 



(30) 



V -P/Vl -P/V2 ••• / 

is a matrix of 2 x 2 blocks such that after the transformation ( p9"D the system 
of equations (^) is rewritten as 

1 



-Z' = BZ 



P 



with matrix 

of block- diagonal form: 

B{p) 



B = P-'AP 



( fin 
B 22 



-P-'P 



In' 



(31) 
(32) 



V 







\ 



Bnn I 



(33) 



Expanding 



P = P 
B = B 



-Pi 



B 1 



:P? 



;B 2 + 



(34) 



and solving equation ([32]) in each power of 1/ p one obtains for Bu 

a 3 

B u = XJ + H + -L u + ^-Lu (35) 

p z p J 



where a and 3 are defined in (|14]) and (|15D and 

is 2 x 2 matrix (g|). 

The 2iV x 2iV matrix system (§) is equivalent to the set of 2 x 2 matrix 
systems of the form 

-Z' = BZ (37) 

P 

where B is one of the Bu. 

Now we will construct solutions of the system of ordinary differential 

equations (f?T|). Let us take 



/' 



p . p-y Y -H --*2 



Then differential equation for U will be 

U' = VU (39) 

where 

V -({a + j I ) h ( a+ i + fi) (Jl + y)+" < 40 ' 

The last system of equations could be transformed to diagonal form because 
the determinant of matrix D Q = T>\ ^ ^ 0. The eigen values of matrix 



V 



1 
ali 



are v\ 2 — ±\/ a h = ±^- Therefore taking the transformation matrix 



r=(\ _4_ 

V 2 2v 



one finds the form of transformed system 

W = TU 
£ = T^VT 

W = SW (41) 



with 



1/ £k + ^(h + k) £k-^(h + k) \ 
+ P {-£k-^(h + k) -£k + ^(h + k) r 



Mi + k) ( \ -\ 

P \ 2 2 



+ ... (42) 



Matrices £i )2 which stand by 1/p and 1/p 2 correspondingly could be trans- 
formed to diagonal form as well after one more substitution 

W = VQ 

Q' = TQ (43) 

Expanding V and T in powers of 1/p one could find the form of matrix JF 
solving matrix equation 

T = V~ X SV - V~ X V (44) 



The result is 



where 



jr = jr + Ijr 1 + ... (45) 

P 



?o=: _°, (46) 



n- ° Zh + spih + h)) m 



-1/ 

2^2 + — 

v £k + ^(h + k) 

Finally we are able to write down the asymptotic solution of system 

Qx =explup+ l^-k^ — (^1 + ^2)) hip> 

Q 2 = expl-up + i-^-k + ^-^1 + k)) Inpj (48) 
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where 



a 



(2N)\ 



l -rf--r 



P 



|^(l+^)((2(l-/i) 2 




(49) 



Using this expression we could restore the asymptotic solution of initial ma- 
trix system of differential equations (|9|) and therefore the asymptotic solution 
of equation ([I]) could be written in the form 

tf(p)=p°$i(p) + p 6 $ 2 (p) (50) 

where $i,2(p) are some one- valued functions and a, b are given by the expres- 
sions 



4^ + ^ (i ' +w y 



(51) 



The coefficients /i and l 2 (see fl28|) and fl36|)) are expressed through the 
elements i^ of the matrix inverse to T (([24|) and (|27|) ). Therefore in order 
to determine the asymptotic expression for solution of (|TJ) we just need to 
calculate the inverse matrix T _1 . 

Matrix T has the form 

/ 1 1 1 ... 1 \ 



T 



Ai 
A? 



A. 



Af + Ai 



A 2 



U 



2AT-1 \2JV-1 



2JV-1 \2N-2 \2AT-1 



; Af 



A TV + 2 

Ai + A 



v2AT-l 



iV 



2AT-2 \2JV-2 



A V H 2~ A W / 



(52) 
Its determinant will not change if we subtract odd columns from even ones: 



det T = det T\ 



1 
Ai 

A? 




l 

2 

Ai 



1 

A 2 



K 



2JV-1 2iV-l\2iV-2 \2iV-l 



^1 




1 

2 
A AT 

2JV-1 \2iV-2 
2 A N 



(53) 
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The last expression could be written in the following form 



detT = ^ d^.-d^W]^ 



(54) 



where W is Wandermond determinant 

1 

W 



1 
A 2 



/'i 

2.V-J ..2N-1 \2JV-1 



1 
UN 



X{' ' m" ' A5" ••• /^ 



27V-1 



(55) 



Using the well known rule of calculation of Wandermond determinant we 
obtain for det T 

Implementing the same trick for the calculation of corresponding minors My 
we obtain the expression for the elements of inverse matrix T -1 : 



f 



1 



1JV 



rife=2(^fe - A 



E 



1 



N 

n 



1 



A i ~~ Al fe=2 ( Afc — Al / 



JV 



t 2 7V - -2 IJ 7T T^ 

fc =2 l A fc _ A L 



(57) 



Then 



/!+/ 2 



JV 1 

k=2 ( A fc _ A l)' 
N I V -1 

-2 ^ — - — n — - — - 

3=2 A i ~ A i fc=2 ( A fc _ A i)' 



(58) 



The product and the sum entering last equation could be evaluated using 
the characteristic equation fl20|) for matrix A ([H]). Indeed, if one write the 
characteristic polynomial P 2 v i n the form 



P2N — (A — A x ). ..(A — A 2 jv) 
which for N large enough become approximately 

P 2 v~(A-A 1 ) 2 ...(A-A iV ) 2 



(59) 
(60) 
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because all the eigen values of matrix Aq become twice degenerate, as it is 
explained earlier, then 



^Ia=a 1 ^(Ai-A 2 ) 2 ...(A 1 -A^ 



(61) 



exactly the required product. Taking into account expression (|^) for P 2 n in 
the limit iV — > oo one easily finds 



A' 



n(Afc-Ai 



(2N)\ (C 2 \ /CAi 
' ' cos < 



fc=2 



(tcy N v 4 



(62) 



In a similar way 



N 



1 



^A \ 
k=2 A fe - Ai 



,k=2 i.i.rn=2,rn^k 



(Al — A m ) 



and then 



N 



rifc=2(Ai - a a 



i 



\ 



p. 



("1) 



2iV 



P, 



(2) 



2/i 



A=Ai 



k=2 Afc — Ai 2y2 



(63) 



(64) 



If one substitute these expressions into (|5T|) one finds that the powers of 
p which describe the branching of solution at infinity become 



a± 



2(1 - a) 2 - M 2 / m 2 I ,, 2n2 2/ 2 

± ^((l-/i 2 ) 2 - M 2 m 2 



C\/2 J(l - /i) 2 - M 2 /m 2 (2V2 



2(1 - a) 2 - M 2 /m 2 If, 0,0 9 / 9 \ , x 

; ± F ((l-/i 2 ) 2 -M7m 2 (65) 



0A/(l-//) 2 -M 2 /m 2 C2V2 



The plus or minus signs are taken at infinities in S + and S- components of 
Sf correspondingly 

3 Mass spectrum 

Having determined the form (p0|) of asymptotic solutions of equation ([!]) at in- 
finities in S + and S_ components of the Riemannian surface Sp we can try to 
impose the condition of decreasing of the wave function at both infinities. But 
the asymptotic expression for the solution fl5"0| ) contains branching functions 
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if powers a and b are not integer. This means that the asymptotic expression 
( |50|) could not be a good approximation for the one valued wave function 
\1/. It is known from complex analysis, that if the powers a and b of p in 
asymptotic expression ([50]) for the wave function \& are not integer numbers, 
the approximation ( [50"D is valid only in some sector (pi < Arg(l/p) < <p 2 
in the neighbourhood of infinity point. Different approximations are valid 
in different overlaping sectors. This behaviour of the asymptotic approxi- 
mation is called Stokes phenomenon ||. But if both a and b are integers, 
the expression (|50"D is a good approximation for non-branching wave function 
everywhere in the punctured neighbourhoods of infinities. 
Due two the double degeneracy of eigen values of matrix A$ the infinite ma- 
trix equation is reduced to the infinite set of second order matrix equations. 
This second order equations can be analysed in the same way as we do this 
in nonrelativistic quantum mechanics. And we conclude therefore that the 
wave function with needed boundary conditions on the both sides of the Rie- 
mannian surface exists only in the situation described above - when a± and 
b± are integers. 

The conditions for a, b to be integers are equivalent to 

2(1 -/i) 2 -M 2 /m 2 

n 



CV2^/(1 -fi) 2 - M 2 /m 2 
^((l-^-M 2 /m 2 )^k (66) 

where n and k are integers. These conditions define the mass spectrum for 
the equation (p. 

The the motion of dust shell in its own gravitational field is parameterised 
by three parameters - the bare mass of the shell M, the Schwarzchild mass 
which is measured at "right" infinity m = m out and the Schwarzchild mass 
inside the shell (or the mass at "left" infinity) m in = fim. There are two 



equations (|66|) on these three parameters, which we could solve for two of 
them. Therefore, the spectrum of Schwarzchild mass m is parameterised by 
two discrete parameters n and k and one continuous (for example /i). But if 
we restrict ourselves with situations when the gravitational field in which the 
shell moves is produced only by the shell itself, that is there is no space-time 
singularity inside the shell, then /i = and the mass spectrum of such states 
is discrete. 
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The expression fl56|) for the mass spectrum is valid only when (1 — p) 2 — 
M 2 /m 2 > when the square root expression entering ( |66D is well defined. 
This restriction corresponds to classical restriction on the trajectories of 
bound motion of the shell. If it is satisfied, then there exist some maxi- 
mal radius of expansion of the shell p(t) < p max . 

Otherwise the shell could escape to infinity in its classical motion. For 
such shells first of conditions (J66|) disappears and the only left is 



^((l-^f-Mym^k (67) 

It is remarkable that for light-like shells (M = 0) in their own gravita- 
tional field (/x = 0) this quantization condition become (we recall that 

2 
pi 

m = m p ;v2k (68) 



C = ml/ 2m 2 



4 Conclusion 

In this technical note we developed a method of constructing asymptotical 
solutions (^) of relativistic Schroedinger equation in finite differences ([[]), 
which describe a self-gravitating spherically symmetric dust shell of matter 
[|1|] . This finite difference equation is defined over Riemannian surface Sp (f|) . 
The configuration space of the system is the real section of this surface and 
has nontrivial geometry. This complicated geometry of configuration space 
of quantum system is due to complicated geometric structure of classical 
space-time manifold of self-gravitating shell. Analysis of the behaviour of 
asymptotic solution on Riemannian surface Sp enabled us to write down 
quantization conditions ( |S6] ) under which the one valued solution of (fl|) with 
physical boundary conditions exists. These quantization conditions define 
discrete spectrum for Schwarzchild mass m if the mass inside the shell is zero. 
The main result of this paper is that this spectrum up to the numerical factors 
coincides with the spectrum, found in [|lj in large black hole approximation 
and in in quasiclassical approximation. This discrete mass spectrum is of 
the form 

m ~ rripivk (69) 

for the light-like shells (shells with bare mass M = 0). This form of quantum 
black hole spectrum was extensively discussed by many authors [|], 0. 
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